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A  Rapid  and  Inexpensive  PCR-Based  STR  Genotyping 
Method  eor  Identieying  Forensic  Specimens 


INTRODUCTION 

The  possibility  exists  that  samples  recovered  from  a  fatal 
accident  site  can  be  misidentified.  Results  from  routine 
toxicological  or  other  biochemical  testing  of  samples  that 
conflict  with  the  preliminary  identiflcation  suggest  this 
possibility  but  are  not  deflnitive.  Genotyping,  a  DNA- 
based  analysis,  which  is  not  a  standard  part  of  sample 
evaluation,  can  serve  as  an  additional  and  independent 
test  for  sample  identiflcation.  Capillary  electrophoresis 
is  the  method  often  used  for  analysis  of  the  polymerase 
chain  reaction  (PCR)  products  generated  in  genotyping 
tests  (1).  This  method  requires  expensive,  dedicated  in¬ 
strumentation,  software  for  analysis,  and  trained  personnel 
potentially  available  only  in  specialized  facilities.  There 
is  a  need  for  a  more  accessible  and  inexpensive  human 
identity  testing  method  in  instances  where  genotyping 
is  necessary  only  as  a  confirmatory  test  and  there  are  a 
minimal  number  of  samples. 

Short  tandem  repeat  (STR)  loci  are  polymorphic 
regions  found  in  the  genome  that  are  used  as  genetic 
markers  for  human  identity  testing.  Typing  of  STR 
loci  by  polymerase  chain  reaction  (PCR)  is  becoming  a 
standard  for  nuclear  DNA  genotyping  analysis.  The  FBI’s 
Combined  DNA  Index  System  (CODIS)  has  defined  a 
core  set  of  STR  markers  for  human  identity  testing  (2) 
which  are  highly  discriminatory  over  a  wide  range  of 
ethnicities.  Robust  primer  sequences  for  these  loci  have 
been  made  publicly  available  (3)  and  are  available  online 
through  the  National  Institute  of  Standards  and  Testing 
(NIST)  Website,  STRbase  (www.cstl.nist.gov/div831/ 
strbase/)  (4).  The  probability  of  samples  from  different 
subjects  having  identical  alleles  for  all  of  the  CODIS  loci 
is  very  small,  leading  to  the  use  of  these  loci  in  forensic 
evaluations  (1,  5). 

Results  will  be  presented  here  describing  a  qualita¬ 
tive  method  of  genotyping  using  the  Agilent  2100 
Bioanalyzer  as  the  microfluidics  platform  for  separation 
of  PCR  products  amplified  from  STR  loci  selected  from 
tetrameric  CODIS  STR  loci.  The  study  demonstrates 
the  novel  use  of  this  instrument  for  rapid,  inexpensive 
electrophoresis  and  analysis  of  STR  PCR  products  to 
allow  discrimination  among  samples  from  aviation  ac¬ 
cidents.  Accuracy  and  reproducibility  were  demonstrated 
with  control  samples  from  a  variety  of  tissues  and  blood. 


The  method  was  applied  to  an  aviation  case  study  where 
initial  sample  identiflcation  did  not  correlate  with  sub¬ 
sequent  biochemical  and  toxicological  evaluation.  A 
paternity  group  sample  set  previously  used  for  capillary 
electrophoresis  (CE)  system  validation  also  was  tested. 
Finally,  the  findings  were  validated  on  a  CE  system  with 
agreement  in  all  cases. 

MATERIALS  AND  METHODS 

Genomic  DNA  extraction.  Tissue  or  blood  was  stored 
frozen  at  -20“C.  Blood  was  thawed  and  a  200  pi  aliquot 
used  for  genomic  DNA  isolation  with  the  QIAamp 
DNA  mini  kit,  #51304  (Qiagen;  Valencia,  CA)  using  the 
protocol  provided  by  the  manufacturer.  Tissue  samples 
included  kidney,  lung,  liver,  muscle,  and  skin.  A  scalpel 
was  used  to  cut  and  mince  approximately  25  mg  from 
the  frozen  tissue.  The  minced  tissue  was  processed  as 
directed  by  the  manufacturer  using  the  ChargeSwitch 
gDNA  MiniTissue  kit,  #CS  11204  (Invitrogen;  Carlsbad, 
CA).  Yield  was  determined  by  A260  on  a  NanoDrop 
Spectrophotometer  (Model  ND- 1 000,  NanoDrop  Tech¬ 
nologies;  Wilmington,  DE).  Quality  was  examined  on  a 
Bioanalyzer  DNA7500  chip  (Agilent Technologies;  Palo 
Alto,  CA).  Samples  showing  presence  of  material  equal 
to  or  greater  than  lOKbp  were  used  for  PCR. 

Sample  sets.  The  1 5  control  samples  were  from  three 
individuals.  SamplesCl-C5  werefromSubject  1,  Samples 
C6-C10  from  Subject  2,  and  Samples  C11-C15  were 
from  Subject  3.  Each  control  sample  set  contained  a 
sample  from  blood,  liver,  lung,  kidney  and  muscle.  The 
aviation  case  study  contained  five  samples,  labeled  TI¬ 
TS.  They  were  from  kidney,  skin,  kidney,  muscle,  and 
lung,  respectively.  The  paternity  set  was  made  up  of  four 
blood  samples,  1304-1307  from  the  College  of  American 
Pathologists,  AABB/CAP  PARF  Survey  B  2005. 

PCR  Primers  and  Amplification.  Sequences  for  the 
ten  primer  pairs  used  in  the  study  were  as  published  (3) 
with  modifications  as  shown  at  STRbase,  www.cstl.nist. 
gov/biotech/strbase/promega_primers.htm  (4).  Nine 
of  these  are  from  the  13  CODIS  core  STR  loci  (6)  and 
one,  amelogenin,  is  a  sex-determination  locus  (5).  Table 
1  lists  the  STR  loci  and  primer  sequences  used  for  am¬ 
plification  and  the  expected  product  sizes.  Primers  were 
synthesized  (Integrated  DNA  Technologies;  Coralville, 
LA)  and  diluted  as  needed. 
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Table  1.  Genetic  loci  and  associated  primer  pair  sequences  used  for  PCR  in  the  study  (4  ).  See 
Materials  &  Methods  for  PCR  conditions. 


Locus* 

Primer  A 

Primer  B 

Product  size 
range  (bp) 

Repeat 

Structure 

Amel 

ACC  TCA  TCC  TGG 
GCA  CCC  TGG  TT 

AGG  CTT  GAG  GCC 
AAC  CAT  CAG 

212,  218 

NA 

CSF1PO 

AAC  CTG  AGT  CTG 
CCA  AGG  ACT  AGC 

TTC  CAC  ACA  CCA 
CTG  GCC  ATC  TTC 

291-327 

[AGATJn 

D3S1358 

ACT  GCA  GTC  CAA 

TCT  GGG  T 

ATG  AAA  TCA  ACA 
GAG  GCT  TGC 

111-143 

TCTA[TCTG]n[T 

CTA]n 

D5S818 

GGT  GAT  TTT  CCT 

CCT 

TGG  TAT  CC 

AGC  CAC  AGT  TTA 
CAA  CAT  TTG  TAT  CT 

123-155 

[AGATJn 

D7S820 

ATG  TTG  GTC  AGG 
CTG  ACT  ATG 

GAT  TCC  ACA  TTT 

ATC  CTC  ATT  GAC 

211-251 

[GATAJn 

D13S317 

ACA  GAA  GTC  TGG 
GAT  GTG  GAG  GA 

GGC  AGC  CCA  AAA 
AGA  CAG  A 

157-201 

[TATCJn 

D16S539 

GGG  GGT  CTA  AGA 
GCT  TGT  AAA  AAG 

GTT  TGT  GTG  TGC 

ATC  TGT  AAG  CAT 
GTATC 

264-304 

[GATAJn 

THOI 

ATT  CAA  AGG  GTA 

TCT  GGG  CTC  TGG 

GTG  GGC  TGA  AAA 
GCT  CCC  GAT  TAT 

171-215 

[AATGJn 

TPOX 

ACT  GGC  ACA  GAA 
CAG  GCA  CTT  AGG 

GGA  GGA  ACT  GGG 
AAC  CAC  ACA  GGT 

TA 

220-256 

[AATGJn 

vWA^ 

GCC  CTA  GTG  GAT 
GAT  AAG  AAT  AAT 

CAG  TAT  GTG 

GGA  CAG  ATG  ATA 
AAT  ACA  TAG  GAT 
GGA  TGG 

123-181 

[TCTA]n[TCTG]n[ 

TCTAJn 

(TCCA[TCTA]n) 

*Amel  (Amelogenin)  is  a  sex-determination  locus,  all  others  listed  are  tetrameric  STR  loci 
^vWA  repeat  sequence  in  ( )  is  not  always  present 
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Amplification  was  performed  using  reaction  conditions 
recommended  in  the  Promega  PowerPlex  2.1  system  (7) 
except  that  the  buffer  was  lx  Amplitaq  gold  reaction 
buffer,  0.1%TritonX-100,  and  0.2mM  each  dNTP  in  a 
25|J.I  final  volume.  Cycling  was  performed  on  a  GeneAmp 
9600  thermal  cycler  (PerkinElmer;  Wellesley,  MA)  using 
the  recommended  cycling  profile  for  40  cycles: 

An  initial  95°  C  11  min,  96°  C  1  min;  then  94°  C 
30sec,  ramp  68  sec  to  60°  C,  hold  for  30  sec,  ramp  50 
sec  to  70°  C,  hold  for  45  sec,  for  10  cycles.  Then  90°  C 
30  sec,  ramp  60  sec  to  60°  C,  hold  for  30  sec,  ramp  50 
sec  to  70°  C,  hold  for  45  sec,  for  30  cycles.  Then  60°  C 
for  30  min,  4°C  soak. 

Preliminary  testing  was  performed  to  optimize  the 
reactions  for  primer  and  template  concentration.  The 
conditions  selected  for  the  study  were  l-5ng  template 
DNA  and  20  pmol  each  primer  for  a  25|J.l  reaction. 
For  multiplexed  reactions,  conditions  were  unchanged 
except  that  the  template  concentrations  were  increased 
two-fold. 

STR  Electrophoresis  and  Analysis.  No  cleanup  of  the 
amplification  products  was  performed.  Detection  of 
a  l|ll  aliquot  from  a  1:4  dilution  in  distilled  water  of 
each  amplification  product  was  performed  on  a  2100 
Bioanalyzer  using  Expert  version  B. 01. 02. SI  136  soft¬ 
ware  and  DNA  500  Labchip  kits  (Agilent  Technologies; 
Palo  Alto,  CA).  Expert  software  was  used  to  overlay 
the  electropherograms  and  assign  PCR  product  length. 
Product  lengths  were  calculated  by  the  software  based 
on  product  migration  times  compared  to  spiked-in  size 
standards.  Capillary  electrophoresis  was  performed  on  a 
310  Genetic  Analyzer  (Applied  Biosystems;  Foster  City, 
CA).  Amplification  and  labeling  of  PCR  products  for 
capillary  electrophoresis  were  done  using  the  PowerPlex 
16  System  (Promega  Corp.;  Madison,  WI). 

RESULTS  AND  DISCUSSION 

The  goal  of  this  study  was  to  develop  a  simple  DNA- 
based  protocol  that  would  allow  an  independent  evalua¬ 
tion  of  accident  site  samples  in  cases  where  toxicological 
testing  results  did  not  agree  with  the  assigned  sample 
identity.  PCR  is  the  recognized  method  of  choice  because 
it  requires  small  amounts  of  starting  DNA  template,  is 
highly  reproducible,  and  is  relatively  insensitive  to  deg¬ 
radation  for  the  amplicon  sizes  used  in  this  study.  An 
extensive  human  genotyping  system  has  been  established 
using  PCR  of  well-characterized,  short  repetitive  genomic 
regions  with  standardized  primers.  The  PCR  primers  were 
selected  from  the  well-characterized  FBI  human  CODIS 
STR  (short  tandem  repeats)  loci.  STRbase  was  used  as  the 
source  of  the  CODIS  primer  sequences  (www.cstl.nist. 


gov/div83I/strbase/)  (4).  An  Agilent  2100  Bioanalyzer 
was  used  for  the  electrophoresis  and  analysis  of  PCR 
samples  as  a  less  expensive,  convenient  alternative  to 
capillary  electrophoresis .  The  Bioanalyzer  uses  microfluid¬ 
ics  chips  to  perform  electrophoresis  in  the  presence  of  a 
fluorescent,  intercalating  dye.  In  the  study  presented  here, 
this  instrument  was  used  for  the  separation  and  analysis 
of  unpurified,  unlabelled  tetrameric  STR  PCR  products 
from  three  sets  of  samples;  first,  a  control  panel  of  tis¬ 
sue  and  blood  samples  from  three  individuals;  second,  a 
case  study  where  toxicological  testing  results  conflicted 
with  sample  identification  made  at  the  accident  site;  and 
finally,  a  paternity  sample  set  that  had  previously  been 
examined  by  capillary  electrophoresis  (CE). 

Nine  primer  sets  from  the  CODIS  STR  loci  set  (6)  and 
agender  identification  locus,  amelogenin  (5)  were  selected 
for  PCR  to  have  an  acceptable  power  of  discrimination. 
Method  development  was  performed  using  the  control 
panel.  The  control  panel  consisted  of  liver,  lung,  kidney, 
muscle,  and  blood  from  each  of  three  subjects.  DNA  was 
extracted  from  all  1 5  samples  and  the  optimum  primer 
and  template  concentrations  were  determined  with  the 
standard  thermocycling  protocol  for  PCR  given  in  the 
Materials  and  Methods  section. 

There  were  two  primary  concerns  with  the  develop¬ 
ment  of  the  assay:  reproducibility  across  tissue  types  at 
each  locus,  and  assay  sensitivity.  Efficacy  of  the  Bioana¬ 
lyzer  for  fragment  sizing  has  been  demonstrated  (8)  with 
DNA7500  chips,  which  have  a  100-7500  basepair  (bp) 
optimum  range.  This  study  showed  reproducible  sizing 
results  for  both  well-to-well  and  day-to-day  replicates 
using  a  single  380  bp  PCR  product  and  a  commercial 
plasmid  digest.  The  DNA  500  chips  used  in  this  study 
have  a  specified  size  range  of 25-500  bp  and  at  least  a  five 
basepair  resolution.  The  CODIS  loci  used  in  the  study  are 
tetrameric.  This  defined  the  resolution  necessary  for  the 
STR  assay  to  be  four  basepairs.  Table  1  lists  the  loci  and 
primers  selected  for  the  study  and  the  expected  product 
size  ranges.  In  the  preliminary  tests,  a  commercial  25  bp 
ladder  containing  1 1  species  was  used,  as  well  as  control 
sample  PCR  products,  to  more  carefully  examine  the 
reproducibility  and  resolution  of  the  DNA  500  chip.  The 
electropherograms  from  the  test  showed  that  there  was 
a  high  level  of  intra-chip  and  inter-chip  reproducibility. 
Table  2A  shows  that  the  CV  (coefficient  of  variation)  of 
the  well-to-well  basepair  values  for  the  25  bp  ladder  were 
low,  ranging  from  0-0.46%  over  6  replicates  of  1 1  spe¬ 
cies/well.  The  size  of  the  50  bp  species  was  called  the  least 
accurately,  showing  a  10%  difference  from  the  expected 
size  (45  bp  instead  of  50  bp).  However,  species  in  the  size 
range  of  the  STR  loci  used  in  this  study  (approximately 
100-300bp)  averaged  approximately  2%  larger  than  the 
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Table  2B.  Comparison  of  the  two  amelogenin  PCR  products  with  sizes  in  basepairs.  %CV 
is  used  to  compare  intrachip  results,  t-test  results  compare  peaks  1  and  2  between  chips. 


Well 

Chip  1 

Peaki 

Peak2 

Peaki 

Chip  2 

Peak2 

2 

205 

209 

206 

210 

3 

205 

209 

206 

210 

4 

204 

209 

206 

210 

5 

205 

209 

206 

210 

6 

205 

209 

207 

210 

7 

204 

208 

207 

210 

8 

204 

208 

206 

210 

9 

203 

208 

206 

210 

%CV 

0.24 

0.22 

0.22 

0.00 

t-Test  (between 
_ chips) _ 

Peak  1  Peak  2 
0.000383  0.000136 


expected  sizes.  The  amelogenin  PCR  products  from  the 
control  panel  also  showed  a  low  CV  of  0-0.24%  for  peaks 
1  and  2  (Table  2B).  This  indicated  a  good  well-to-well 
correlation.  Paired  t-Tests  were  performed  to  compare 
data  between  chips  and  showed  statistically  significant 
differences  for  both  the  25  bp  ladder  and  the  amelogenin 
PCR  products  (p<0.05,  Tables  2A  and  B). 

Accuracy  of  the  Bioanalyzer  depended  on  the  size  of  the 
molecule  but  was  well  within  specifications  for  the  chip. 
The  25bp  ladder  species  were  called  consistently  within 
2%  of  the  expected  size  but  were  generally  called  larger 
than  the  expected  size  (Table  2A),  except  for  the  50bp 
species,  which  was  consistently  called  smaller.  Table  2B 
shows  that  the  amelogenin  PCR  products  were  consis¬ 
tently  called  within  4%  of  the  expected  size  by  the  Expert 
software.  These  have  been  published  as  being  212  and 
218  bp,  respectively  (5)  but  migrate  consistently  as  seven 
to  nine  basepairs  smaller  in  our  hands  on  the  Bioanalyzer 
platform.  The  difference  in  migration  may  be  because  the 
PCR  reaction  was  not  purified  or  denatured  before  elec¬ 
trophoresis  or  to  the  presence  of  “ragged”  single  stranded 
ends  on  the  PCR  products,  compared  to  the  25bp  ladder 
species  in  the  commercial  product.  This  suggested  that 
the  Bioanalyzer  was  limited  in  the  detection  of  absolute 
sizes  of  the  genotyping  PCR  products.  However,  the 
high  level  of  reproducibility  of  the  intra-chip  size  calls 
also  suggested  that  reliable  comparative  analysis  would 
be  possible  and  would  allow  an  accurate  determination 
of  forensic  sample  identification. 

Reproducibility  between  different  samples  from  the 
same  subject  was  evaluated.  DNA  isolated  from  control 


samples  C3  and  C5  (isolated  from  lung  and  muscle)  was 
subjected  to  PCR,  and  the  products  were  analyzed  on 
the  Bioanalyzer  in  replicate  wells.  No  significant  differ¬ 
ence  in  fragment  lengths  called  by  the  Expert  software 
was  seen  when  the  replicates  were  compared  across  tissue 
types  (data  not  shown).  This  suggested  that  the  isolation 
protocol  yielded  consistent  template  material  regardless 
of  starting  tissue. 

Based  on  our  sample  identification  need,  the  design  of 
the  assay  was  for  rapid  evaluation  and  relative  comparison 
of  PCRproducts,  rather  than  absolute  identification  ofthe 
specific  alleles  displayed  by  the  samples.  Due  to  statistical 
significance  of  inter-chip  sensitivity,  all  loci  comparisons 
were  done  between  samples  electrophoresed  on  the  same 
chip  (Table  2). 

All  1 5  of  the  control  samples  were  used  for  PCR  with 
the  ten  primer  pairs  using  the  conditions  listed  in  Materi¬ 
als  and  Methods  and  subjected  to  electrophoresis  on  the 
Bioanalyzer.  The  results  were  evaluated  using  the  Expert 
software.  The  electrophoresis  results,  both  aligned  sizes 
(in  basepairs)  and  electropherograms,  were  examined  for 
patterns  to  discriminate  between  samples.  It  was  noted 
that  the  occurrence  of  singlets  (homozygous  alleles)  versus 
doublets  (heterozygous  alleles)  was  sufficient  to  discrimi¬ 
nate  between  some  samples  such  that  the  three  control 
subjects  (samples  C1-C15)  could  be  identified  based  on 
the  pattern  of  homozygous  alleles  alone.  This  is  a  useful 
starting  point  for  data  analysis;  however,  homozygosity 
was  not  necessary  to  determine  sample  identities  since 
differences  at  biallelic  loci  were  observed. 
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The  aviation  case  study  contained  five  samples,  T 1 -T5 . 
PCRwas  performed  using  the  established  conditions  with 
each  of  the  ten  primer  sets,  followed  by  electrophoresis 
on  the  Bioanalyzer.  The  results  showed  thatT5  had  ho¬ 
mozygous  alleles  for  the  CSFIPO,  D5S818,  andTHOl 
loci.  This  was  sufficient  to  identify  Sample  T5  as  being 
unique  fromTl  -T4  (Table  3A) .  The  basepair  sizes  assigned 
to  the  PCR  products  suggested  that  the  remaining  four 
samples  could  be  divided  into  two  groups.  Samples  T 1 
and  T3  fell  into  one  group  and  Samples  T2  and  T4  fell 
into  a  second.  Samples  in  each  group  potentially  belong 
to  the  same  individual  based  on  the  assigned  basepair  sizes 
for  the  PCR  products  (Table  3A).  Examining  overlays  of 
the  electropherograms  confirmed  that  SamplesTl  andT3 
displayed  the  same  peaks  for  all  loci  tested  and  differed 
fromT2  andT4.  Samples  T2  andT4  also  displayed  the 
same  peaks  for  all  loci  tested.  This  confirmed  the  assign¬ 
ment  to  three  groups.  Figure  1  is  an  electropherogram 
overlay  for  the  vWA  locus  PCR  products  from  samples 
T1-T5,  showing  an  example  of  the  relative  peak  migra¬ 
tion  for  the  aviation  case  study  samples.  Figure  1 A  clearly 
shows  the  difference  in  migration  pattern  between  peaks 
from  samples  Tl,  T2,  andT5.  Figure  IB  shows  electro¬ 
phoretic  migration  patterns  clearly  demonstrating  that 
samples  Tl  andT3  contained  the  same  vWA  alleles  that 
differed  from  samples  T2  andT4,  which  in  turn,  display 
the  same  vWA  alleles.  SeeTable  3A  for  the  corresponding 
assigned  basepair  measurements  and  assessment.  Com¬ 
bining  the  human  identity  testing  information  with  the 
toxicological  data  suggested  that  the  original  assignment 
ofsamplesT2,T4,  andT5  to  the  same  individual  was  an 
error.  Sample  T5  was  clearly  from  a  different  individual. 
Samples  Tl  and  T3  were  blind  controls  from  the  same 
individual  included  to  confirm  the  efficacy  of  the  assay 
for  quality  assurance. 

To  further  validate  the  assay,  the  aviation  case  study 
samples  were  examined  by  CE  (Materials  and  Methods) 
to  confirm  the  Bioanalyzer  conclusions,  and  a  set  of  four 
samples  making  up  a  paternity  group  (samples  1304- 
1307),  which  had  previously  been  examined  by  CE,  was 
tested  using  the  Bioanalyzer.  The  CE  results  are  sum¬ 
marized  in  Table  3,  along  with  the  Bioanalyzer  basepair 
assignments  and  analysis  summary.  Listed  in  Table  3A  are 
the  results  of  the  aviation  case  study  by  sample  with  both 
the  base  pair  size  assigned  by  the  Expert  software  (“bp” 
column)  and  the  allele  assignment  made  by  the  Geno typer 
software  following  capillary  electrophoresis  on  an  ABI 
310  (“al”  column).  Results  were  summarized  using  only 
the  microfluidics  chip  electrophoresis  data  and  agreed 
with  the  CE  results  at  all  loci  except  CSF 1 PO  and  Tho  1 . 
Alleles  11  and  12  for  locus  CSFIPO  were  not  resolved 
sufficiently  to  be  called  with  the  Bioanalyzer  software. 
Manual  examination  of  the  electropherograms  showed 


a  partially  resolved  doublet.  In  the  studies  reported  here, 
this  did  not  change  the  final  conclusion  for  the  samples. 
The  poor  resolution  of  these  alleles  is  thought  to  be  due 
to  the  larger  sizes  of  the  CSFIPO  PCR  products.  Alleles 
11  and  12  are  reported  as  311  and  315  bp  respectively. 
Ffeterozygous  samples  from  all  other  loci  used  in  the 
study  with  alleles  differing  by  four  basepairs  were  easily 
resolved.  For  Bioanalyzer  genotyping,  it  is  expected  that 
choosing  loci  with  PCR  products  of  less  than  300bp  will 
have  better  sensitivity  and  accuracy.  Alleles  9  and  9.3  for 
the  Tho  1  locus  could  not  be  resolved.  Even  with  manual 
observation,  a  doublet  call  could  not  be  made.  There  is 
only  a  three  basepair  size  difference  between  these,  195 
and  198  bp,  respectively.  However,  as  shown  in  this 
study,  it  is  unlikely  that  an  inability  to  resolve  these  two 
alleles  will  affect  the  overall  accuracy  of  the  assay  across 
a  wide  range  of  loci. 

The  paternity  group  samples  were  known  to  be  a 
mother  (1304),  daughter  (1305),  and  potential  fathers 
(1306  and  1307).  The  microfluidics  chip  results  (Table 
3B)  clearly  identified  four  individuals,  and  did  not  sup¬ 
port  1307  as  a  potential  father  but  did  support  paternity 
for  1306  for  all  loci  tested.  These  findings  agreed  with 
the  CE  results. 

The  ability  of  the  Bioanalyzer  to  resolve  products 
from  multiplexed  PCR  reactions  was  tested.  The  limited 
resolution  of  the  Bioanalyzer  chip  and  unlabeled  PCR 
products  prevented  the  high  degree  of  multiplexing 
possible  with  the  CE  systems.  However,  a  series  of  five 
duplex  reactions  were  successfully  resolved.  This  reduced 
reagent  and  microchip  usage  by  half  CSFIPO  andThoI 
reactions  could  be  cycled  and  electrophoresed  together,  as 
could  D13S17  andTPOX,  D3S1358  and  amelogenin, 
D16S539  and  D5S818,  and  finally,  vWA  and  D7S820 
(see  Material  and  Methods  for  conditions)  with  no  loss 
of  resolution  and  no  overlap  in  the  electrophoresed  PCR 
products.  Figure  2  shows  an  example  of  this  duplex  reac¬ 
tion  set.  Secondary  products  can  be  seen  in  the  D5S818 
and  D16S539  panel.  However,  the  height  and  migration 
time  clearly  show  these,  as  well  as  the  shoulder  seen  in 
the  TPOX  peak,  as  artifacts.  The  duplexed  loci  were 
combined  for  optimum  separation  but  are  not  the  only 
multiplex  reactions.  A  number  of  additional  combina¬ 
tions  were  tested  (data  not  shown),  and  there  appeared 
to  be  no  excessive  secondary  products  or  inhibition  of 
amplification  due  to  the  presence  of  multiple  sets  of 
primers  at  the  concentrations  reported  here.  Indeed,  a 
triplex  of  D5S818,  CSFIPO,  andTHOl  amplifies  and 
separates  well. 

To  address  the  sensitivity  of  the  assay  to  contaminated 
template  samples  that  could  result  in  a  miscalled  geno¬ 
type,  PCR  products  from  a  mixture  of  two  templates 
were  examined.  Samples  1304  and  1306  were  mixed 
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Table  3.  Bioanalyzer  (bp)  and  capillary  electrophoresis  (al)  results  for  samples  in  the  Case  Study  and  Paternity  Group.  The  Bioanalyzer  summary 
was  supported  by  the  CE  results  in  all  cases,  see  text. 


B.  The  paternity  group  results  with  summary  and  paternity 
determination  from  the  Bioanalyzer  results. 
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A.  The  aviation  case  study  results  with  summary  determined  from  the 
Bioanalyzer  results. 
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Table  3  (continued) 


*bp  basepair  from  the  Bioanalyzer  assignment,  al  allele  from  the  CE  results 
tNA  not  applicable 
JND  not  discriminatory 


Figure  1.  Electropherogram  overlays  of  locus  vWA  PCR  products.  The  X-axIs 
is  migration  time,  the  Y-axis  is  fluorescence  units. 

Overlay  of  samples  T 1 ,  T2,  T5  show  their  unique  migration  patterns. 

Overlay  of  samples  T1,  T3  and  T2,  T4  respectively,  show  matching  migration 
patterns,  suggesting  T1=T3  and  T2=T4. 
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Figure  2.  Duplex  PCR  electropherograms  for  the  ten  STR  loci  used  in  the  study  for  Control 
sample  2.  PCR  performed  using  conditions  described  in  Materials  and  Methods.  Peaks  with 
intersecting  horizontal  lines  were  considered  valid  products  for  analysis  by  the  Bioanalyzer 
software.  The  X-axis  is  migration  time,  the  Y-axis  is  fluorescence  units. 
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in  three  proportions  (Fig.  3A).  Both  the  mixed  and  un¬ 
mixed  samples  were  subjected  to  PCR  with  D3S1358 
primers  and  the  products  subjected  to  electrophoresis 
on  the  microfluidics  chip,  as  above.  Control  samples  Cl 
and  C14  were  similarly  mixed  and  PCR  was  performed 
with  D7S820  primers.  In  Figure  3  the  electrophoresis 
profiles  from  both  of  these  series  are  shown.  Figure  3A 
shows  the  presence  of  additional  peaks  or  shoulders  in 
the  1304/1306  mixtures,  clear  evidence  of  more  than 
one  template  source.  However,  the  mixture  of  samples 
Cl  and  C14  shown  in  Figure  3B  has  a  confusing  profile 
that  could  result  in  a  false  genotype  assignment  in  both 
the  75:25  and  50:50  C14:C1  mixtures  due  to  the  more 
intense  homozygous  Cl  peak.  The  significantly  lower 
left-most  peak  contributed  by  C14  and  the  shoulder  on 
the  Cl  peak  (also  contributed  by  Cl 4)  would  suggest 
the  possibility  of  contaminated  template.  Therefore, 
care  should  be  taken  to  look  for  greatly  differing  peak 
heights  as  indicators  of  sample  contamination.  These 
results  suggest  that  only  heavily  contaminated  samples 
will  result  in  confusing  results  since  contaminating  peaks 
could  not  be  seen  in  10:90  and  90:10  mixtures  of  these 
samples  (data  not  shown).  Taken  together.  Figures  3A 
and  B  suggest  that  in  the  event  of  contaminated  samples 
an  ambiguous  reading  would  likely  be  clearly  diagnosed 
by  results  at  other  loci. 

The  statistical  power  of  the  assay  was  estimated  us¬ 
ing  data  from  the  Einum  and  Scarpetta  study  for  three 
North  American  ethnic  groups  (9).  The  nine  STR  loci 


were  determined  to  provide  a  combined  matching 
probability  of  approximately  9.5E+07  to  3E+08  if  all 
tetrameric  alleles  could  be  detected  with  the  Bioana¬ 
lyzer  protocol.  Removing  the  CSFIPO  locus  from  the 
matching  probability  calculation  naturally  reduces  the 
matching  probability  (Table  4)  but  the  assay  retains  an 
acceptable  level  of  power  for  small,  defined  populations. 
The  assay  was  not  able  to  discriminate  between  the  fairly 
common  9  and  9.3  Thol  alleles  (see  Table  3).  Binning 
of  these  two  alleles,  using  the  formula  for  power  of 
discrimination  (10),  provided  a  correction  for  the  Thol 
locus  (M.  Scarpetta,  personal  communication)  giving  a 
corrected  matching  probability  of  approximately  5E+07 
to  1.9E+08  (Table  4).  It  has  been  observed  that  the 
Thol  9.3  and  10  alleles  that  differ  by  one  basepair  are 
difficult  to  resolve  on  some  platforms  (11,  12).  Allele 
10  is  relatively  rare  and  has  not  been  detected  to  date 
in  our  samples  sets,  but  a  binning  of  alleles  9.3  and  10 
(performed  as  for  the  9  and  9.3  allele  binning)  showed 
a  small  change  in  allelic  frequencies  from  0.8372  with 
no  binning  to  0.8355  for  the  Black  population,  0.8873 
to  0.8824  for  the  Caucasian  population,  and  0.8295 
to  0.8272  for  the  Hispanic  population.  Binning  these 
alleles  would  have  only  a  minimal  effect  on  the  overall 
matching  probability  of  this  assay.  The  amelogenin 
gender  discrimination  locus  also  was  included  in  the 
study  and  does  provide  additional  power,  but  because  it 
is  not  an  STR  locus,  it  was  not  included  in  the  Einum 
and  Scarpetta  data. 


Table  4.  Estimation  of  Matching  Probability  for  the  nine  STR  loci  used  in  the 
study  based  on  PD  values  from  Einum  and  Scarpetta  (9). 


Matching  Probability  (1/1-PD)^ 

STR  locus 

N.  American  Blacks 
N=8000 

N.  American  Caucasians 
N=8000 

N.  American 
Hispanic  N=1000 

CSF1PO* 

12.66 

8.26 

8.06 

D5S818 

9.71 

7.14 

8.26 

D7S820 

12.05 

16.13 

13.51 

D13S317 

7.35 

13.33 

17.86 

D16S539 

14.49 

11.36 

12.82 

THOfl 

9.71 

11.90 

8.62 

TPOX 

11.63 

5.32 

6.76 

vWA 

17.24 

10.31 

12.99 

Total  MP 

3.07E+08 

9.42E+07 

1.56E+08 

w/o  CSF1 PO 

2.43E+07 

1.14E+07 

1.93E+07 

w/oTHOI  &  CSF1PO 

2.50E+06 

9.57E+05 

2.24E+06 

‘Cannot  discriminate  all  tetrameric  alleles  for  this  locus 
^PD  Power  of  Discrimination,  see  text 

*  Includes  correction  for  inability  to  discriminate  alleles  9  and  9.3,  see  text. 
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Figure  3.  The  appearance  of  PCR  products  from  mixture  of  two  sample  templates  in  three  proportions. 

A.  Appearance  of  D3S1358  PCR  products  from  mixed  template  Paternity  group  samples  1304  and  1306  in 
three  proportions  where  the  presence  of  both  templates  is  easily  detected.  Peaks  1  and  3  are  from  Sample 
1304  and  Peaks  2  and  4  are  from  Sample  1306,  as  indicated  in  the  overlay  panel. 

B.  Appearance  of  mixed  samples  following  PCR  with  D7S820  using  Control  samples  Cl  and  C14.  The 
overlay  shows  the  relative  mobility  of  the  peaks  from  each  sample  with  peaks  1  and  2  from  Cl  4  and  peak  3 
from  Cl .  See  text  for  discussion. 
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